Surface polaritons on left-handed spheres 
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We consider the interaction of an electromagnetic field with a left-handed sphere, i.e., with a 
sphere fabricated from a left-handed material, in the framework of complex angular momentum 
techniques. We emphasize more particularly, from a semiclassical point of view, the resonant aspects 
of the problem linked to the existence of surface polaritons. We prove that the long-lived resonant 
modes can be classified into distinct families, each family being generated by one surface polariton 
propagating close to the sphere surface and we physically describe all the surface polaritons by 
providing, for each one, its dispersion relation and its damping. This can be achieved by noting 
that each surface polariton corresponds to a particular Regge pole of the electric part (TM) or 
the magnetic part (TE) of the S matrix of the sphere. Moreover, for both polarizations, we find 
that there exists a particular surface polariton which corresponds, in the large radius limit, to that 
supported by the plane interface. There also exists, for both polarizations, an infinite family of 
surface polaritons of whispering gallery type having no analogs in the plane interface case and 
specific to left-handed materials. They present a "left-handed behavior" (phase and group velocities 
are opposite) as well as a very weak damping. They could be very useful in the context of plasmonics 
or cavity quantum electrodynamics. 

PACS numbers: 78.20.Ci, 41.20.Jb, 73.20.Mf, 42.25. Fx 
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I. INTRODUCTION 

In the present paper, we continue our analysis of the 
scattering of electromagnetic waves by objects of sim- 
ple shape made of a left-handed material. (We refer 
to Refs. 1.2,3.4.5.6 for important articles dealing with 
such materials.) Our analysis is based on the complex 
angular momentum (CAM) method which allows us to 
completely describe the resonant aspects of the problem. 
In a previous work, we have studied scattering from a 
left-handed cylinder (disk)*. In the present paper, we 
shall extend this analysis to scattering from a left-handed 
sphere because of its numerous potential practical appli- 
cations. 

Scattering from left-handed spheres has been already 
considered in some recent article o 8 ' 9 : 10 : 11 ' 12 : 13 : 14 : 15 : 16 ' 17 
and the aspects linked to the existence of the resonant 
surface polariton modes (RSPM's) supported by left- 
handed spheres have been more particularly considered 
in Refs. UHlJ. In the present paper, by using CAM 
techniques in connection with asymptotics beyond all 
order o 18 : 19 : 20 ' 21 , we shall provide a clear physical explana- 
tion for the excitation mechanism of the RSPM's of the 
sphere as well as a simple mathematical description of 
the surface waves, i.e., of the so-called surface polaritons 
(SP's), that generate them. We refer to the Introduction 
of Ref. [22J for a presentation of the CAM method and 
for a short bibliography as well as to the monographs of 
Newton 2 ^, Nussenzveig2* and Grandy** for more details. 
It should be noted that all the techniques which we use 
here are well-known techniques of electromagnetism of 
ordinary dielectric media. It is only recently that they 



have been introduced in the context of electromagnetism 
of dispersive media (see Refs.Ml22.26). They enable us to 
go well beyond the work completed until now to describe 
the SP's propagating on objects of simple shape. 

Our paper is organized as follows. Section II is de- 
voted to the exact theory: we introduce our notations, 
we provide the expression of the S matrix of the system 
and we then discuss the resonant aspects of the prob- 
lem for both polarizations. In Sec. Ill, by using CAM 
techniques, we qualitatively describe the SP's supported 
by the left-handed sphere and we establish the connec- 
tion between these SP's and the associated RSPM's. In 
Sec. IV, by using asymptotic techniques, we describe 
semiclassically the different SP's and we provide analytic 
expressions for their dispersion relations and their damp- 
ing. We show more particularly the existence of SP's of 
whispering gallery type. Finally, in Sec. V, we conclude 
our paper by emphasizing the main results of our work 
and by briefly discussing the implication of some of our 
results in the context of plasmonics and cavity quantum 
electrodynamics. 



II. EXACT S MATRICES AND SCATTERING 
RESONANCES 

A. General theory 

From now on, we consider the interaction of an elec- 
tromagnetic field with a sphere of radius a fabricated 
from a metamaterial and having an effective frequency- 
dependent permittivity e(ui) and an effective frequency- 



2 



dependent permeability (J-(lo). Here, and in the following, 
we implicitly assume the time dependence exp(— iu/t) for 
electric and magnetic fields. We consider that the sphere 
is embedded in a host medium of infinite extent having 
the electromagnetic properties of the vacuum. We intro- 
duce the usual spherical coordinate system (r, 8, lo). It is 
chosen so that the sphere and surrounding medium re- 
spectively occupy the regions corresponding to the range 
< r < a (region II) and to the range r > a (region I). 
Furthermore, in order to describe wave propagation, we 
also introduce the wave number 



c 



(1) 



where c denotes the velocity of light in vacuum, and the 
refractive index of the sphere 



n(u>) = \J e(w)n(uj). 



(2) 



As far as the electric permittivity e(u>) and the mag- 
netic permeability fi(co) of the sphere are concerned, we 
choose the expressions already used in Ref.0for the cylin- 
der. We assume that they are respectively given by 



and 



e(uo) 



Flo 2 



(1-F) 



(3) 



(4) 



where < F < 1 and uob = luq/v 1 — F. We are aware 
that these expressions do not describe the electromag- 
netic behavior of all the real left-handed media but they 
are the most used in the literature. Moreover, even if we 
considered more complicated expressions for the electric 
permittivity and magnetic permeability, our main results 
would remain valid (see the Conclusion of Ref. 0). 

For the theoretical aspects of our work, we shall as- 
sume that luq < uji, < uj p . We have e(w) < in the 
frequency range u) € ]0, ui p [ and /u(oj) < in the fre- 
quency range u> S ]wo,Wb[. Thus, the electric permittiv- 
ity, the magnetic permeability and the refractive index 
are simultaneously negative in the region loq < lo < 
In that region, the metamaterial presents a left-handed 
behavior. As far as the numerical aspects of our work 
are concerned, we shall work with F = 0.4 and with 
the reduced frequencies Lo^afc = 5.52, Lo^a/c m 7.127 and 
u> p a/c = 11.04. Even though we restrict ourselves to that 
configuration, the results we shall obtain numerically are, 
in fact, very general and they permit us to correctly il- 
lustrate the theory. Furthermore, we have chosen these 
particular values, which have been already used in Ref. 7 
for the cylinder, in order to be able to compare the three- 
dimensional problem with the two-dimensional one. 

The S matrix of the sphere is of fundamental impor- 
tance because it contains all the information about the 
interaction of the sphere with the electromagnetic field. 



It can be obtained from Maxwell's equations and usual 
continuity conditions for the electric and magnetic fields 
at the interface between regions I and I I 24 i 25 ' 27 . Because 
of the spherical symmetry of the scatterer, the S ma- 
trix is diagonal and its elements Sw(lo) are given by 
See>(uj) = Si(lo) Sw. For our problem, the elements of 
the electric part (TM polarization) of the S matrix are 
given by 



with 



Sf(u) = l-2af(u) 



af (to) 



Df(u)' 



(5) 



(6) 



where Cf{w) and Df(to) are two 2x2 determinants 
which are explicitly given by 



Cf (u) — \J e (lo) I (wj'ipt [n (u>) u>a/c\ ij)' e (wa/c) 

—ipi(ioa/c)ip'i[n(Lo)ioa/c\ (7a) 

Df(uj) = \J e (lo) (u>)ip£ [n (to) toa/c] c}^ (ioa/c) 

— Q 1 (wa/e)^[n(w)wa/c], (7b) 

while the elements of its magnetic part (TE polarization) 
are given by 



with 



Sf(w) = l-2af(w) 



(8) 



(9) 



where Cf 1 (lo) and Df 1 (lo) are also two 2x2 determinants 
which are explicitly given by 



Cf 1 (lo) — \J e (lo) I fi (to)ipe (locl/c) [n (lo) Loa/c] 

~ipi[n(Lo)Loa/c\ip[(Loa/c) (10a) 

Df(Lo) = yje (lo) In (w)d 1} (coa/c) $ [n (lo) Loa/c] 



~ipi[n(Lo)Loa/c\C t ^ (uia/c) 



(10b) 



In Eqs. dTj) and (flTJ)) . we use the Ricatti-Bessel func- 
tions ijje(z) and C^\z) which are linked to the spherical 
Bessel functions je(z) and h^(z) by ipi(z) — zjz(z) and 

<W(z) = zh ( e 1) (z) (see Ref.il). 

From the S matrix elements, we can, in particular, 
construct the total scattering cross section of the sphere 
which is defined as the ratio of the total energy scattered 
and the incident energy intercepted. It can be expressed 
in terms of the coefficients af(uj) and af 1 (lo) and it is 
given by^^l 



<jt(lo) 



2tt 



[fcM] 



^£(2^ + l)[|af(a;)| 2 + K M (a;)| 2 ]. (11) 



From the S matrix elements, we can also precisely de- 
scribe the resonant behavior of the sphere as well as the 



3 



geometrical and diffractive aspects of the scattering pro- 
cess. Here, the dual structure of the S matrix plays a cru- 
cial role. Indeed, the S matrix is a function of both the 
frequency u> and the angular momentum index £. It can 
be analytically extended into the complex u> plane as well 
as into the complex A plane (the CAM plane). (Here, A 
denotes the complex angular momentum index replacing 
£+1/2 where £ is the ordinary momentum inde x 23 i 24 ' 25 .) 
The poles of the S matrix lying in the fourth quadrant 
of the complex u) plane are the complex frequencies of 
the resonant modes. These resonances are determined 
by solving 



Df' M (u) = for 



e n*. 



The solutions of (fT2")l are denoted by 



,(°) 



(12) 
iTi p /2 

where > and Ti p > 0, the index p permitting us to 
distinguish between the different roots of Q12p for a given 
In the immediate neighborhood of the resonance u>£ P , 
(uj) has the Breit-Wigner form, i.e., proportional to 

Tip/2 



yE.M 



UJ — UJ 



(°) 



iT lp /2 



(13) 



The structure of the S matrix in the complex A plane 
allows us, by using integration contour deformations, 
Cauchy's theorem and asymptotic analysis, to provide 
a semiclassical description of scattering in terms of rays 
(geometrical and diffracted). In that context, the poles 
of the S'-matrix lying in the CAM plane (the so-called 
Regge poles) are associated with diffraction. They are 
determined by solving 



D A-1/ 2 M 



for uj > 0. 



(14) 



Of course, when a connection between these two faces of 
the S'-matrix can be established, the resonance aspects 
of scattering are then semiclassically interpreted. 



B. Numerical results and physical interpretation 

In Fig. QJi, we display the total cross section plotted 
as a function of the reduced frequency uiajc. Rapid vari- 
ations of sharp characteristic shapes can be observed. 
Such a strongly fluctuating behavior is due to the scat- 
tering resonances associated with the long-lived resonant 
modes of the sphere, i.e., the long-lived resonant states 
of the photon-sphere system: when a pole of the S ma- 
trix is sufficiently close to the real axis in the complex 
uj plane, it has a strong influence on the total cross sec- 
tion. In Fig. [T]d, resonances are exhibited for the two 
polarizations. For certain frequencies, we can clearly ob- 
serve a one-to-one correspondence between the peaks of 
<jt(uj) and the resonances near the real ua/c axis but, 
in general, the situation seems very confused. This is 
due to the profusion of long-lived resonant modes in and 
around the frequency range where the sphere presents a 
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FIG. 1: a) Total cross section or. b) Scattering resonances 
in the complex uia/c plane. Dots (•) correspond to poles of 
S e (uj) while plus (+) correspond to poles of S M (lo). 



left-handed behavior. In spite of that, it should be nev- 
ertheless noted that (i) the strongly fluctuating behav- 
ior generated by electric resonances (TM polarization) is 
localized within and slightly around the frequency range 
where the sphere presents a left-handed behavior and (ii) 
the strongly fluctuating behavior generated by magnetic 
resonances (TE polarization) is totally localized within 
that frequency range. Furthermore, by zooming in on 
the distribution of resonances in regions close to the real 
axis of the complex u> plane, we have also observed accu- 
mulations of resonances for large values of I: 

(i) For the TM polarization, there exists an accu- 
mulation of resonances which converges to the limiting 
frequency uj s satisfying 



e{uj s ) + 1 = 



and given by 



(15) 



(16) 



We have for the corresponding numerical reduced fre- 
quency uj s a/c ~ 7.806. 

(ii) For the TE polarization, there exists an ac- 
cumulation of resonances at the limiting frequency ujf 
satisfying 



H(cj f ) + 1 = 



and given by 



UJf = UJ 



2-F 



(17) 



(18) 



We have for the corresponding numerical reduced fre- 
quency ujfd/c « 6.172. 
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(iii) For both polarizations, there exists an accu- 
mulation of resonances at the pole loq of /i(lo) which cor- 
responds more precisely to 



III. SEMICLASSICAL ANALYSIS: REGGE 
POLES, SURFACE POLARITONS, AND 
RESONANCES 



-00. 



(19) 



From now on, we shall more particularly focus our at- 
tention on the physical interpretation of the long-lived 
resonant modes whose excitation frequencies belong to 
frequency ranges in which n(u>) < 0. We shall prove 
that the corresponding resonances are generated by ex- 
ponentially small attenuated SP's propagating close to 
the sphere surface and we shall provide a numerical and a 
theoretical description of these surface waves. The corre- 
sponding resonant modes are the so-called RSPM's. For 
the other resonant modes associated with bulk polari- 
tons (the resonant modes whose excitation frequencies 
belong to the frequency range in which n(u>) > and 
those with shorter lifetime), we do not provide a similar 
analysis. This is not very serious as they do not have, in 
physics, the importance of RSPM's. Indeed, in the new 
field of plasmonics, those are the SP's with long propaga- 
tion lengths and therefore very small attenuations that 
are especially interesting from the point of view of the 
practical applications. Furthermore, if we consider the 
photon-sphere system as an artificial atom for which the 
photon plays the usual role of the electron, we must then 
keep in mind that, in the scattering of a photon with fre- 
quency l>£° , a decaying state (i.e., a quasibound state) of 
the photon-sphere system is formed. It has a finite life- 
time proportional to l/Te p . The resonant states whose 
complex frequencies belong to the family generated by 
SP's are therefore the most interesting because they are 
very long-lived states. 

We shall finally conclude this section by making a 
brief comparison with the results obtained in our pre- 
vious studies concerning the left-handed cylinder— and 
the metallic and semiconducting spheres 25 : 

(i) The total cross sections for the left-handed 
cylinder and the left-handed sphere as well as their spec- 
tra of resonances are rather similar. However, it should 
be noted that in the scattering by a sphere both the TM 
and TE polarizations contribute to the total cross sec- 
tion (we recall that for the cylinder, the two polarizations 
can be studied separately). But, it is also important to 
note that, for both scatterers, SP's and their associated 
RSPM's correspond to only one polarization (the TE po- 
larization for the cylinder and the TM polarization for 
the sphere). 

(ii) As far as the spectrum of resonances is con- 
cerned, the left-handed sphere is a physical system much 
richer than the metallic or the semiconducting sphere (see 
Figs. 1 and 2 of Ref. |H and the discussion at the end 
of Sec. II of that reference) and this is certainly very 
interesting for practical applications of left-handed clcc- 
tromagnetism and more particularly in plasmonics and 
in cavity quantum electrodynamics. 



A. Complex angular momentum approach: Results 

In the CAM approac h 23 i 24 ' 25 , Regge poles determined 
by solving Eq. (TT4")) are crucial to describe diffraction as 
well as resonance phenomena in terms of surface waves. 
From the Regge trajectories associated with the SP's sup- 
ported by the left-handed sphere, i.e., from the curves 
Asp = Asp(w) traced out in the CAM plane by the corre- 
sponding Regge poles as a function of the frequency, we 
can more particularly deduce the following: 

(i) the dispersion relations of SP's which connect 
their wave numbers ksp with the frequency to: 



fcsp(w) 



Re Asp(w) 



(20) 



(ii) the dampings ImAsp(w) of SP's, 

(iii) the phase velocities v p as well as the group ve- 



locities v g of SP's: 



P Re Asp (w) 



and v„ 



d Re Asp(w) ' 



(21) 



(iv) the semiclassical formula (a Bohr-Sommerfeld- 
type quantization condition) which provides the location 

of the excitation frequencies wSl, of the resonances gen- 
erated by SP's: 



Re A 



sp 



( W $p) =1+1/2, I: 



1,2,. 



(22) 



(v) the semiclassical formula which provides the 
widths of these resonances 



esp 



Im Asp (lo) [d/duj Re Asp (w)] 



[d/dujRe A S p(w)] 2 + [d/du Im A S p(w)] 2 



(23) 

and which reduces, in the frequency range where the con- 
dition \d/dwKe Asp(w)| ^S> \d/dujlm Asp(w)| is satisfied, 
to 



Im Asp(w) 



d/du Re Asp (w) 



(24) 



All these results can be established by generalizing, mu- 
tatis mutandis, our approach and our calculations devel- 
oped in Refs. |7||22| for dispersive cylinders. The transition 
from two dimensions to three dimensions induces some 
additional technical difficulties (vectorial treatment, ex- 
istence of a caustic, asymptotics for spherical harmonics, 
etc.) which can be overcome following and extending the 
works of Newton in quantum mechanics (see Chap. 13 of 
Ref. |23T ) and the works of Nussenzveig^ and Grandy 25 
in electromagnetism of ordinary dielectric media. 
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FIG. 2: Regge poles in the CAM plane for the TM polariza- 
tion, a) The distribution corresponds to uoa/c — 5.6. b) The 
distribution corresponds to coa/c = 6.2. c) The distribution 
corresponds to uoa/c — 7.0. 



FIG. 3: Regge poles in the CAM plane for the TE polariza- 
tion, a) The distribution corresponds to uoa/c — 5.6. b) The 
distribution corresponds to uoa/c = 6.2. c) The distribution 
corresponds to uoa/c — 7.0. 



B. Regge poles and Regge trajectories 

Figures [5] and [3] exhibit the distribution of Regge poles 
for both polarizations for three different reduced frequen- 
cies lying in the frequency region where n(uo) < 0. We 
have identified and indicated some particular Regge poles 
which are associated with surface waves orbiting around 
the left-handed sphere and which explain its resonant be- 
havior. These figures are at first sight qualitatively sim- 
ilar to Figs. (4) and (5) of Ref. H where we displayed the 
corresponding Regge pole structures for the left-handed 



cylinder. It should be however noted that, in the tran- 
sition from the cylinder to the sphere, the TM and TE 
polarizations are exchanged and curvature corrections in- 
duce modifications in the quantitative behavior (as we 
shall show in Sec. IV). More precisely, for both polariza- 
tions, we can observe in Figs. [2] and [3] the following. 

- One of these Regge poles is associated with the 
SP noted SPoo which, as we shall show in Sec. IV, cor- 
responds in the large radius limit (i.e., for a — > oo) to 
a SP supported b y the pla ne interface and theoretically 
described in Refs.lUlSlM 
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FIG. 4: Regge trajectory for the Regge pole associated with 
SPoo (TM polarization). 



FIG. 6: Regge trajectory for the Regge pole associated with 
SPoo (TE polarization). 
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FIG. 5: Regge trajectories for the Regge poles associated with 
the first three whispering-gallery SP's (TM polarization). 



FIG. 7: Regge trajectories for the Regge poles associated with 
the first three whispering-gallery SP's (TE polarization). 



- The other Regge poles are associated with an 
infinite family of SP's of whispering gallery type denoted 
by WGSP„ with n <E W which have no analogs for the 
plane interface (see Sec. IV) as well as for curved metallic 
or semiconducting interfaces (see Refs. I22ll26l ). 

In Figs. HE we have displayed the Regge trajectories 
of the first SP's for the TM and TE polarizations. They 
have been obtained by solving numerically Eq. (fT4"|) . We 
can observe some interesting features: 

- The dispersion curve for the surface wave SPoo 
of the TM polarization is a positive and monotonically in- 



creasing function of u>. As a consequence, the associated 
group and phase velocities given by Eq. (f2Tj) are both pos- 
itive and SPoo has an ordinary behavior. It should be also 
noted that this SP exists in the frequency range ui G]0, u> s [ 
and therefore in the range ]ljo,^&[ where the refraction 
index is negative but also outside this range. For low val- 
ues of oj, its damping becomes very large and thus this 
surface wave has a negligible role for these frequencies in 
the scattering process and in the resonance mechanism. 
Furthermore, as u> — » u> s the dispersion curve increases 
indefinitely. This result will permit us to explain the ac- 
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FIG. 8: Resonances generated by SP^ (TM polarization). 



FIG. 10: Resonances generated by SP^ (TE polarization). 
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FIG. 9: Resonances generated by the first three whispering 
gallery SP's (TM polarization). 



FIG. 11: Resonances generated by the first three whispering- 
gallery SP's (TE polarization). 



cumulation of resonances which converge to the limiting 
frequency uj s for the TM polarization. 

- The dispersion curve for the surface wave SPoo 
of the TE polarization is a positive and monotonically 
decreasing function of u>. As a consequence, the associ- 
ated phase velocity is positive while the group velocity 
is negative [see Eq. (pH]) ]. SPoo has a "left-handed be- 
havior" . It should be also noted that this SP only ex- 
ists in the frequency range ]ujf,uj}y[ which is included in 
the frequency range ]wo, Wf,[ where the refraction index 
is negative. Its damping is always weak and thus this 
surface wave always plays a significant role in the scat- 
tering process and in the resonance mechanism. Finally, 
as w — > ujf the dispersion curve increases indefinitely. 
This result will permit us to explain the accumulation of 
resonances which converge to the limiting frequency Wf 
for the TE polarization. 

- As far as the surface waves WGSP„ (n e N*) 
of the TM and TE polarizations are concerned, they 
present, at first sight, a behavior which is rather inde- 
pendent of the polarization. The real part of a given 
Regge pole Awgsp„ vanishes for a frequency in the range 



]u>o, u>b[ and becomes negative. This Regge pole then mi- 
grates to the third quadrant of the CAM plane and be- 
comes unphysical. So we can consider that the surface 
waves WGSP„ with n 6 N* only exist in a subdomain of 
the frequency range ]luq, lui,[ where the refraction index is 
negative. The dispersion relations of all these surface 
waves are positive but monotonically decreasing func- 
tions. Their group velocities are always negative while 
their phase velocities are positive [see Eq. (|2ip]. All these 
SP's thus have a left-handed behavior. Furthermore, be- 
cause the dampings of these surface waves are always 
weak, they all play a significant role in the scattering pro- 
cess and in the resonance mechanism. Finally, it should 
be noted that as u> — ► ojq, the dispersion curves increase 
indefinitely. This result will permit us to explain the ac- 
cumulation of resonances which converges to the limiting 
frequency luq for the TM and TE polarizations. 
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FIG. 12: Zoom-in on the total cross section ot- 



C. Resonances 

In Figs.l8l fTT1 we present samples of complex frequencies 
for the RSPM's associated with the surface waves SPoo 
and WGSP„ with n = 1, 2 and 3. They have been cal- 
culated from the semiclassical formulas (|2"2"|) and (f2"3")l by 
using the Regge trajectories determined numerically by 
solving Eq. (TTi|) (see Figs. |3J[7]). A comparison between 
the semiclassical spectra and the "exact ones" (calculated 
by solving numerically Eq. 112ft shows a very good agree- 
ment. Moreover, we can also observe some interesting 
features (mutatis mutandis they were already present for 
the left-handed cylinder): 

- The resonance spectrum associated with the 
surface wave SPoo of the TM polarization (see Fig. [8|) 
extends beyond the frequency range ]wo,W(,[ where the 
sphere presents a left-handed behavior because SPoo ex- 
ists for oj s]0, w s [. Furthermore, inserted into the semi- 
classical formulas ([22j) and (|23|) . the behavior of the 
Regge trajectory of Asp^ near oj s easily explains the ex- 
istence of the family of resonances close to the real axis 
of the complex u) plane which converges for large £ to the 
limiting frequency ui s . 

- The resonance spectrum associated with the sur- 
face wave SPoo of the TE polarization (see Fig. nU|) fully 
lies inside the frequency range ]ct>o,w&[ where the sphere 
presents a left-handed behavior because SPoo exists only 
in that range. Furthermore, inserted into the semiclas- 
sical formulas (|2"2"|) and (|23p . the behavior of the Regge 
trajectory of Asp^ near u)t explains the existence of the 
family of resonances close to the real axis of the com- 
plex u> plane which converges for large I to the limiting 
frequency lu f . 

- The resonance spectra associated with the sur- 
face waves WGSP„ (n e N*) of the TM and TE polar- 
izations (see Figs. l9l and [TTj) fully lie inside the frequency 
range ]ujq, u)b[ where the sphere presents a left-handed be- 
havior because all these surface waves exist only in that 
range. Furthermore, inserted into the semiclassical for- 



mulas (|22p and (|23[) . the behavior of the Regge trajectory 
of a given Regge pole Awgsp„ near ujq explains the ex- 
istence of a corresponding family of resonances close to 
the real axis of the complex uj plane which converges for 
large £ to the limiting frequency ujq. 



In conclusion, we have established a connection be- 
tween the complex frequencies of the long-lived resonant 
modes (or RSPM's) of the left-handed sphere and the 
SP's noted SPoo and WGSP„ with n G N* which are 
supported by its surface. In other words, in spite of the 
great confusion which seems to prevail in the resonance 
spectrum of the left-handed sphere (see Sec. II), we have 
been able to fully classify and physically interpret the 
resonances thanks to CAM techniques. We now invite 
the reader to look at Fig. [T2] where we have zoomed in 
on Fig. [TJ On the total cross sections ctt(w) we have 
identified the peaks corresponding to resonances and, for 
each one, we have specified the SP which has generated 
it as well as its polarization and the associated "quantum 
number" I. This has been achieved by using the results 
displayed on Figs. IBTTTTl It is also possible to identify 
the "quantum numbers" of the different resonant modes 
by plotting the associated distribution of electromagnetic 
energy density. Indeed, the number of local maxima is 
twice the quantum number I. These distributions can 
been obtained from the theory developed in Refs. [321I331 . 
For example, in Figs. IT5HTH1 we have displayed, for the 
TE polarization, the energy distributions for the reso- 
nant modes £ = 9 and £ = 10 generated by SPoo and 
for the resonant modes £ — 8 and £ — 9 generated by 
WGSPi. For obvious symmetry reasons, we have only 
considered one quarter of the equatorial section of the 
sphere. 
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FIG. 13: Electromagnetic energy density for the resonant 
mode I = 9 generated by SPoo (TE polarization). 



FIG. 15: Electromagnetic energy density for the resonant 
mode I = 8 generated by WGSPi (TE polarization). 



SP^ M0 TE 
coa/c = 6.2022 




WGSP, 1=0 TE 




FIG. 14: Electromagnetic energy density for the resonant 
mode £ = 10 generated by SPoo (TE polarization). 



FIG. 16: Electromagnetic energy density for the resonant 
mode 1 = 9 generated by WGSPi (TE polarization). 



IV. ASYMPTOTICS FOR SURFACE 
POLARITONS AND PHYSICAL DESCRIPTION 

In order to obtain a deeper physical understanding of 
the SP's orbiting around the left-handed sphere and to 
justify the terminology previously used (i.e., the nota- 
tions SPoo and WGSP n for the SP's), we must "analyti- 
cally" solve Eq. (fT4| for A = Agp or equivalently 



(lu) Cisp-i/2 <W C ) _ ^A S p-i/2 [n(uj)auj/c] 



V Vi") Cist-1/2 ( auJ / c ) ^Asp-i/2 [n(ui)au/c] 
for the TM polarization and 



(25) 



H{w) Clsp-i/2 ( auJ / c ) _ ^Asp-i/2 [n(u>)aui/c] 
e ( w ) Ci 1 s ) p _i/2 (auj/c) ^Asp-i/2 [n(uj)auj/c] 



(26) 



for the TE polarization. In other words, we seek to obtain 
explicit formulas for the Regge poles of the problem. We 
have previously done it for metallic and semiconducting 
object o 22 i 26 as well as for left-handed cylinders- and we 



shall here extend such an analysis to left-handed spheres. 
We were recently aware of a beautiful article by Berry 34 
where this problem has also been considered for a curved 
dielectric interface with an electric permittivity negative 
and independent of u>. 

By replacing Ricatti-Bessel functions by spherical 
Bessel functions and then using their relations with the 
ordinary Bessel functions (see Ref. HI) 



and 



1/2 

(27) 

as well as Jx(-z) = e mX J x {z) and J' x (-z) = -e t7rX J' x (z) 
(see Ref. \2&j) in order to take into account that n(ui) < 
in the most interesting frequency range, it is easy to prove 
that Eqs. (|2"5"j) and (f2"6")l respectively reduce to 



J'\ SP [\ n (u)\ua/c} 
J\ SP [\n(oj)\ uia/c] 



1 



2n(u 
(uja/c) 



(-) 



+ 



MM h{ 1] {uoa/c) 2 V/ i ( w ) 



(28) 
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for the TM polarization and 



1 



J\ sp [\n(ui)\uia/c] 2n(uj) \u>aJ 



for the TE polarization. These two last equations must 
be compared respectively with Eqs. (37) and (38) of 
Ref. which provide SP Regge poles for the left-handed 
cylinder. The first term on the left-hand side and on the 
right-hand side of Eqs. (|28[) and ([29]) are exactly those 
appearing in Eqs. (37) and (38) of Ref. 0. The remain- 
ing terms are curvature corrections due to the change of 
dimension. Equations Q28p and ([2^]) can be solved follow- 
ing the method used in order to solve Eqs. (37) and (38) 
of Ref. 0, i.e., by using asymptotic analysis. As we have 
previously noted in Ref. |7J, the choice of the asymptotic 
expansions for the Bessel functions strongly depends on 
the relative positions of the arguments \n(u>)\u>a/c and 
uja/c with respect to the complex order Asp. In order to 
simplify the discussion, we choose to describe theoreti- 
cally SP's in the frequency ranges where they generate 
the RSPM's with the longest lifetime (such modes are 
the most important from the physical point of view). In 
other words, we shall seek Asp^ (to) for the TM polariza- 
tion with lu in the neighborhood of lu s , Asp 00 (w) for the 
TE polarization with lu in the neighborhood of LUf and 
Awgsp„ ( w ) for both polarizations with lu in the neigh- 
borhood of qjq . In fact, in spite of these restrictions, we 
shall obtain asymptotic results valid in large frequency 
ranges. 



Asymptotics for surface polaritons of SP C 
and physical description 



-type 



Let us first consider the Regge pole associated with 
SPoo for the TM polarization. We assume lu in the neigh- 
borhood of lo s [but also that u> < lui, so that n(uS) < 
0] and then we can also assume that ReAsp 00 (w) > 



Loa/c > \n(uj)\uja/c and formally that ReAsp oc (w) 3> 1 
and ReAsp DO (w) 3> ImAsp 00 (w). As a consequence, we 
can use the Debye asymptotic expansions for J\(z) and 



(z) valid for large orders (see our previous paper— for 
details as well as Appendix A of Ref. |35| or Ref. \3t 
can then write 



We 



^>(4Wj . [Alp, - (\n{m)\tua/cf 
•h SP ^[\n(uj)\uja/c\ 



,1/2 



\n(Lo)\(uia/ c) 



(30) 



and 



H 



I 1/2 



(coa/c) 



where 



x (i_ie 2Q(Asp »' uo/c) ) (31) 
a(A, z) = (A 2 - z 2 ) 1 ' 2 - A In ( A + (A ' ~ ^ 



(32) 



In Eq. (|3ip , we have taken into account an exponentially 
small contribution (the term exp[2a(Asp oc , Lua/c)]) which 
lies beyond all orders in perturbation theory. This term 
can be captured by carefully taking into account Stokes 
phenomenon and is necessary to extract the asymptotic 
expression of the imaginary part of Agp- [In Eq. (|3ip 
we have given to the Stokes multiplier function the value 
1/2.] For more precision, we refer to our previous paper- 
as well as to Refs.[llll20|l2l|. Now, by inserting © and 
(f3Tj) into Eq. (|28|) . we obtain an equation which can be 
"easily" solved. We first neglect the exponentially small 
term and this equation reduces to a fourth-order poly- 
nomial equation. We only retain the physical solution 
which correctly modifies the formula (47a) obtained in 
Ref. for the corresponding SP of the left-handed cylin- 
der. It provides the real part of Asp 00 (w). We then take 
into account the exponentially small term and we obtain 
perturbatively the imaginary part of Asp^ (lu) . We finally 
have 



ReA SPoo H~ (-) x 



Im Asp^(w) 



lu) - 1 4[1 - \e{cu)\] 2 \uja 



|e(w)[[|e(w)| + K^)] , e 2 H + l / c \2 e 2 (u;)[\e(D)\fj,(u) + l] / c \3 e 2 (uj) re 



' { 1 Hm/ [e 2 (w) -l][l-[e(w)|] 2 \ua) + 4[1 - |e(^)|] 4 W/ 



(-) 



M A „ (Re ^SPooM) 2 - (wtt/c) 2 



(lu) — 1 / Re Asp^ (u 



x exp[2a(Re Asp,^ (w), uia/c) 



(33a) 
(33b) 



Let us now consider the Regge pole associated with we must solve Eq. (|29[) which only differs from Eq. (|28|) 
SPoo of the TE polarization. In order to describe it, by the exchange of e(u>) and h(lu). Here, we assume lu 
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in the neighborhood of uif and then we can also assume 
that Re Asp^ (w) > \n(u))\u>a/c > wa/c and formally that 
ReAsp^M > 1 and ReAsp TO (w) > ImAsp^w). As 
a consequence, we can use again the Debye asymptotic 

expansions for J\(z) and H^\z) valid for large orders 



and the resolution of Eq. (|29|) can be modeled on that 
of Eq. pg]). Formulas and (T5TJ) remain valid and by 
inserting them into Eq. (|29|). we obtain 



ReAsp^^) ~ ( — ) x 



\! m 2 H-i + 4[i-|/xM 



|A*(w)|[|A*(w)| + c(w)] , M 2 H + 1 / c \a ,i'( w )[| M ( w )|e( w ) + l] / c \ 2 , ^(w) 



[At 2 (w) - 1][1 - |a*H|] 2 Vwa/ 4[1 - |m(^)|] 4 \m 



+ 



Im Asp^ (lj) 



M 2 M \ „ (ReA SPoo H) 2 - W<=) 2 



ReA SPoo (w) 



x exp[2a(ReAsp 00 (w),6tia/c)]. 



(34a) 
(34b) 
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FIG. 17: Regge trajectory for the Regge pole associated FIG. 18: Regge trajectory for the Regge pole associated 
with SPoo (TM polarization). Comparison between exact and with SPoo (TE polarization). Comparison between exact and 
asymptotic theories. asymptotic theories. 



Equations (|33aj) - (|33bj) and (|34a|) - (j34b[) provide respec- 
tively analytic expressions for the dispersion relation and 
the damping of the surface polaritons SPoo of the TM 
and TE polarizations. The following important physical 
features must be noted: 

- These two SP's present exponentially small at- 
tenuations. 

- Their wave numbers fcsPoo (w) can be obtained 
from (|33ap and (fSTTj) for the TM polarization and from 
(J34aj) and ([20]) for the TE polarization and they could 
permit us to derive analytically their phase velocity 
u p = w/fcsp 00 (w) as well as their group velocity v g = 
dio/dksp^ (uj). 



- For a — > oo - i.e., in the flat interface limit - the 
wave number ksp^ (w) for the TM polarization reduces 
to 




(35) 



This expression is the usual dispersion relation found in 
Refs. 01(3131] for the p-polarized SP - i.e., the SP for 
which the magnetic field H is normal to the incidence 
plane - supported by the flat interface. For the TM 
polarization, SPoo is therefore the counterpart of the p- 
polarized SP supported by the flat interface. 

- For a — > oo - i.e., in the flat interface limit - the 
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wave number fcgp^ (u) for the TE polarization reduces to 




mMI[ImMI + £H] 

M 2 M - 1 



(36) 



This expression is the usual dispersion relation found in 
Refs. HaHMMI for the s-polarized SP - i.e., the SP for 
which the electric field E is normal to the incidence plane 
- supported by the flat interface. For the TE polarization, 
SPqo is therefore the counterpart of the s-polarized SP 
supported by the flat interface. 

- Furthermore, for a — > oo, the imaginary parts 
(|33b[) and (|34bp of Asp^ vanish. In the flat interface 
limit, SPqo for both polarizations have no damping like 
the SP's supported by the flat interface^^^i. 

- A connection with the results obtained in our 
previous study concerning the left-handed cylinder— can 
be established. By comparing (|33ap with (47a) of Ref. 
and (|34ap with (49a) of Ref. [3 we can notice that SP^ of 
the sphere for the TM (TE) polarization is the counter- 
part of SPoo of the cylinder for the TE (TM) polariza- 
tion. In Eqs. (|33a[) and (|34a[) . the supplementary terms 
are associated with the curvature corrections appearing 
in Eqs. (j2"5|) and (f2"9")l . In certain frequency ranges, these 
corrections are of the same magnitude than the remaining 
terms. As a consequence, it is not possible to obtain dis- 
persion relations on the sphere perturbatively from those 
on the cylinder. 

- The function Re Asp^ (w) given by (|33a[) has a 
pole when e(c<j) + 1 = 0- i.e., for to = lu s . Further- 
more, the imaginary part (|33b[) of Asp^ vanishes for 
lu = lu s . These two results explain semiclassically, for the 
TM polarization, the accumulation of resonances which 
converge to the limiting frequency u> s . 

- The function Re Asp^ (lu) given by (|34a|) has a 
pole when /i(w) + 1 = 0- i.e., for u> = tut. Further- 
more, the imaginary part (|34b[) of Asp^ vanishes for 
lu = u)f. These two results explain semiclassically, for the 
TE polarization, the accumulation of resonances which 
converge to the limiting frequency lu f . 

- We have numerically tested formulas (|33ap and 
((33b)) (see Fig. \T7} as well as formulas (pHaf and p4b| 
(see Fig. [TS]). They provide rather good approximations 
i) for Re Asp^ (lu) in large frequency ranges for both po- 
larizations, ii) for Im Asp^ (lu) in the neighborhood of lu s 
for the TM polarization and iii) for ImAsp oc (w) in the 
neighborhood of LUf for the TE polarization. 



B. Asymptotics for surface polaritons of WGSP„ 
type and physical description 

Let us finally consider the Regge poles associated with 
the surface waves WGSP„ for the TM and TE polar- 
izations. We must now solve Eqs. (|28|) and ([29]) for 
Asp = Awgsp„ by assuming u> in the neighborhood of luq. 
The configuration is identical to that of Ref. |7| for the left- 
handed cylinder. We can consider that Re Awgsp„( w ) > 



Lua/c and formally that Re Awgsp„ ( w ) 3> 1 an d 
ReAwGSP„(^) > ImAwGSP„(^)- Furthermore, A\vgsp„ 
is in the immediate neighborhood of \n(uj)\uja/c i.e., lies 
in the Airy circle centered on \n{us)\u)a/ c. As a conse- 
quence, we can still replace (z) by its Debye asymp- 
totic expansion valid for large orders and therefore insert 
(l3"Tj) into (f2"5)) and (|2"9"j). As far as J\{z) is concerned, 
we must now use its uniform asymptotic expansion (see 
Appendix A of Ref. [HI or Ref. [3a ) 



J x (z)~(2/z) 1/3 Ai \(2/z) 1/3 (A — z) 



(37) 



where Ai{z) denotes the Airy function^. Thus, we must 
insert 



J 



AwGSP„ 



[|n(w)|ci;a/c] 



•AwgspJK w )IWc] 



{2/[\n(co)\Loa/c]} 1/3 



Ai 


{2/[|nM|^a/ C ]} 1/3 [AwGSP„ 


— |n(w)|wo/c] 


Ai 


1 /3 

{2/[\n(u))\ua/c}} [A W gsp„ 


— \n(uj)\uja/c} 



(38) 



into (|28f and (|29|) . We obtain two equations which can 
be solved by using the method already considered in our 
previous work 7 - and which is an adaptation of the method 
invented by Rayleigh a long time ago in order to describe 
mathematically the whispering-gallery phenomenon in 
acoustic o 37 ' 38 (see also Ref. |39|). After having noticed 
that the curvature corrections can be neglected in that 
configuration, we recover the formulas obtained for the 
left-handed cylinder. We have 



Re A 



WGSP 



, (lu) ~ \n(cu)\Lua/c - 



\n(Lu)\Lua j c\ 



(39a) 



Im A 



WGSP., 



e(lu) f \u(lu) | Lua/c\ 2 ^ 3 

1/2 



[Re A W gsp„M] 2 - (Lua/c) 2 



(Lua/c) x n 
x exp [2a(Re Awgsp„( w )> luq/c)] 



(39b) 



for the TM polarization and 



ReAwGSP„(^) ~ \n(Lu)\Lua / c + 



\n(ui)\Lua/ c\ 1 ' 3 



(40a) 



Im A W gsp„ M 



[i(lu) / \u(lu)\ cua/c\ 2 ^ 3 



[Re A W gsp„(w)] 2 - (u>a/cy 



1/2 



(Lua/c) x n 
x exp [2a(Re Awgsp„ ( w )> luq/c)] 



(40b) 
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FIG. 19: Regge trajectory for the Regge pole associated with 
WGSPi (TM polarization). Comparison between exact and 
asymptotic theories. 



der with the polarizations exchanged: (|39a[) and (|39b[) 
must be compared with (55a) and (55b) of Ref. [D while 
(|40a| and (|40bj) must be compared with (55a) and (55c) 
of Ref. 7. Therefore, mutatis mutandis, the physical re- 
marks and analysis of Ref. [7J concerning the whispering- 
gallery surface waves propagating on the cylinder can be 
taken again. We just note that: 

- These SP's have no counterparts in the plane 
interface case. 

- They do not exist for curved metallic or semi- 
conducting interfaces (see Refs. I22ll26l ). 

- The function ReA\yGSP„(w) given by (|39ap or 
(|40ap has a pole which is that of n(w) and therefore 
which corresponds to w — u>q. Furthermore, the imagi- 
nary parts (|39bj) and (|40bj) of Awgsp„ vanish for to = ujq. 
These results justify all our previous remarks concerning 
the accumulations of resonances which converge to the 
limiting frequency ujq. 

- We have numerically tested formulas (|39ap - (|40b() 
(see Figs.[T9landr20lfor WGSPi). They provide very good 
approximations for Re AwGSPi (w) in the full frequency 
range ]u;o,Wf>[ where the sphere presents a left-handed 
behavior. They also provide very good approximations 
for Im AwGSPi (w) in a rather large frequency range above 
the limiting frequency ui . 
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FIG. 20: Regge trajectory for the Regge pole associated with 
WGSPi (TE polarization). Comparison between exact and 
asymptotic theories. 



for the TE polarization. In Eqs. (|39a|) - (|40b[) . we have 
introduced the zeros x n of the Airy function (the first 
three ones are X\ ~ —2.3381..., x 2 ~ —4.0879... and xz ~ 
-5.5205...). 

Equations (|39a[) - (|40b|) provide analytic expressions for 
the dispersion relation and the damping of the surface 
polaritons WGSP„ for both polarizations. These expres- 
sions are those already obtained for the left-handed cylin- 



V. CONCLUSION AND PERSPECTIVES 

In the present paper, we have considered the inter- 
action of electromagnetic waves with a sphere fabricated 
from a left-handed material. We have mainly emphasized 
the resonant aspects of the problem. We have shown that 
the long-lived resonant modes can be classified into dis- 
tinct families, each family being generated by one SP and 
we have physically described all the SP's orbiting around 
the sphere by providing, for each one, a numerical and a 
semiclassical description of its dispersion relation and its 
damping. 

We have more particularly shown that the left-handed 
spherical interface can support both TE and TM polar- 
ized SP's. For each polarization, there exists a partic- 
ular SP which corresponds, in the large radius limit, to 
the SP which is supported by the plane interface and 
which has been theoretically described in Refs. [29ll30ll3ll . 
However, there also exists, for each polarization, an infi- 
nite family of SP's of whispering gallery type and these 
have no analogues in the plane interface case. They 
are only supported by the curved left-handed interface 
(they do not exist for curved metallic or semiconducting 
interface a 22 ' 26 ) and they are analogs to those described 
for the cylindrical interface in Ref. @. 

One could believe at first sight that the existence of 
surface waves of whispering gallery type (and of the as- 
sociated resonant modes) in the frequency region where 
both permittivity and permeability of the sphere are neg- 
ative is rather trivial because they are also present for 
an ordinary dielectric sphere. So, one could also believe 
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that it is inappropriate to consider such surface waves 
as SP's. In fact, that is not at all the case. Indeed, it 

is well known that the excitation frequencies uj^ of the 
whispering-gallery resonant modes of an ordinary dielec- 
tric sphere having a positive and frequency-independent 
refractive index n are obtained from the relation (see e.g. 
Refs. EHH) 

(I + l/2\ 1/3 
n ua/c = £+1/2- ( J x n 

+corrections depending on the polarizationj(41) 

By replacing n by n(w), this relation permits us to de- 
scribe the scattering resonances of the left-handed sphere 
when both permittivity and permeability are positive and 
therefore the scattering resonances appearing on Fig. [1] 
for to > uj p . This result can be obtained by extending the 
approach developed in Refs. [40ll4ll On the other hand, 
one can easily verify that this relation does not provide 
the excitation frequencies of the whispering-gallery res- 
onant modes when the sphere presents the left-handed 
behavior (for wj < w< W(, on Fig. [T]). In this case, it is 
necessary to use the theory which we developed in Sec. 
IV. From (|39a[) or (|40a[) which have been obtained by 
assuming e(w) < and /i(w) < and from (f2"2")) we then 
obtain 

(I + 1/2 \ 1/3 
-n(u)ua/c = £+1/2- I — -^—\ x n + . . (42) 

Equations (l4"T]) and (|4"2"]) seem similar in form because 
they both describe whispering-gallery resonant modes 
but, in fact, they are very different. They have distinct 
ranges of validity and thus describe physical phenomena 
of different nature (right-handed behavior in the first case 
and left-handed behavior as well as accumulation of res- 
onances in the second one). The existence of the left- 
handed whispering-gallery resonant modes is related to 
the conditions e(ui) < and h(oj) < and therefore, 
more particularly (but not only), to the existence of free 



electrons. As a consequence, we think that these modes 
can be called resonant surface polariton modes and that 
the associated surface waves described by (p9]l and (|40|) 
are of SP type. 

Because of the proliferation of SP's, left-handed 
spheres are systems much richer than metallic or 
semiconducting spheres 2 ^. They are therefore much 
more interesting as artificial atoms ("plasmonic 
atoms" ) 22 i 26 ' 42 i 43 and this could have important conse- 
quences in terms of practical applications in the fields of 
plasmonics and nanotechnologies. Here, we shall briefly 
more particularly focus our discussion to the possible 
applications to cavity quantum electrodynamics. It is 
well known that high quality factor whispering-gallery 
modes in microdisks and microspheres, because they are 
associated with a strong localization of the electromag- 
netic field and a very long lifetime for the photons, make 
these resonators interesting for practical applications 
(see, for example, Ref. HU and references therein). They 
have been used, in particular, to produce lasers and 
to enhance the spontaneous emission of quantum dots. 
Up to now, these resonators have been fabricated from 
ordinary solid dielectric materials. We think that the use 
of left-handed material could constitute an interesting 
appealing alternative: indeed, whispering-gallery modes 
supported by left-handed disks or spheres seem to 
present the interesting usual properties (strong localiza- 
tion and very long lifetime) but also very unusual ones 
linked with the left-handed behavior (negative group 
velocities and positive phase velocities). However, at 
this stage of the reflection, it is not possible for us to say 
if quantum electrodynamics in "left-handed resonators" 
could generate new interesting physical phenomena and 
lead to advances in physics. 

Finally, it should be noted that our approach as well 
as our main results are not limited to the left-handed 
materials described by the effective electric permittivity 
and the effective magnetic permeability ([3]) and ([4]) but 
still remain valid for more general left-handed materials 
(see the conclusion of Ref. 0) . 
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